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Spacecraft  Constrained  Maneuver  Planning  for  Moving  Debris 
Avoidance  Using  Positively  Invariant  Constraint  Admissible  Sets 

Avishai  Weiss,  Morgan  Baldwin,  Christopher  Petersen,  R.  Scott  Erwin,  and  Ilya  Kolmanovsky 


Abstract —  This  paper  extends  our  prior  work  on 
spacecraft  debris  avoidance  based  on  positively  invariant 
sets  to  the  case  of  moving  debris.  In  this  approach,  a  time- 
varying  connectivity  graph  is  constructed  between  a  set 
of  forced  equilibria,  forming  a  virtual  net  that  is  centered 
around  a  nominal  orbital  position  of  the  spacecraft.  The 
connectivity  between  two  equilibria  is  determined  based 
on  contractive  sets  in  order  to  guarantee  that  collision- 
free  transitions  between  equilibria  can  be  effected  within 
a  giving  time  period.  A  graph  search  algorithm  is  imple¬ 
mented  to  find  the  shortest  path  around  the  debris,  or 
through  the  debris’  path  once  it  has  moved  on.  Simulation 
results  are  presented  that  illustrate  this  approach. 

I.  Introduction 

Given  the  rapidly  growing  amount  of  debris  in 
Earth  orbit,  spacecraft  maneuvers  must  satisfy  debris 
avoidance  requirements.  While  obstacle  avoidance  has 
been  studied  extensively  in  robotics  applications  [3], 
[4],  the  related  spacecraft  problems  have  several  unique 
features  and  benefit  from  the  development  of  spe¬ 
cialized  approaches  and  algorithms.  In  particular,  the 
space  environment  is  relatively  uncluttered.  Spacecraft 
dynamics  are  different  from  those  of  typical  robots. 
Maneuver  efficiency  with  respect  to  time  and  fuel 
consumption  is  important.  The  states  of  the  spacecraft 
and  the  debris  can  only  be  estimated,  often  with  a 
significant  estimation  error.  Finally,  computational  al¬ 
gorithms  must  be  suitable  for  implementation  on-board 
where  computing  power  is  very  limited. 

In  [1],  we  have  introduced  an  on-board  maneuver 
planning  approach  based  on  the  use  of  constraint- 
admissible  positively  invariant  sets  to  determine  con¬ 
nectivity  between  a  set  of  forced  and  unforced  space- 
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craft  equilibria  forming  a  virtual  net  in  the  vicinity  of 
the  spacecraft.  Two  equilibria  are  connected  if  a  choice 
of  a  Linear  Quadratic  (LQ)  feedback  gain  can  be  made 
that  results  in  a  transition  between  the  equilibria  and 
which  avoids  the  debris  collision  while  satisfying  the 
imposed  thrust  limits.  The  connectivity  graph  for  all  the 
equilibria  in  the  net  is  constructed  and  graph  search 
algorithms  are  employed  to  optimize  an  equilibria 
hopping  sequence  that  avoids  the  debris  collisions.  To 
facilitate  the  on-board  construction  of  the  connectivity 
graph,  a  fast  growth  distance  computation  procedure 
between  two  ellipsoidal  sets  has  been  proposed  in  [1]. 

Unlike  the  open-loop  trajectory  optimization  ap¬ 
proaches,  we  do  not  rely  on  precise  assignment  of 
spacecraft  position  to  the  time  instants  along  the  trajec¬ 
tory,  but  instead  switch  to  the  next  set-point  and  con¬ 
troller  gain  once  appropriate  conditions  are  satisfied. 
While  this  approach  can  be  conservative,  it  facilitates 
fault-tolerant  and  disturbance-tolerant  execution  of  the 
maneuvers.  Furthermore,  by  employing  disturbance- 
invariant  sets  [14],  we  can  assure  robustness  to  unmea¬ 
sured  (but  set-bounded)  disturbances  and  uncertainties. 
This  extension  to  handling  unmeasured  disturbances 
and  uncertainties  using  techniques  of  [14]  is  not  pur¬ 
sued  in  the  present  paper. 

In  [2],  we  incorporated  limited  thrust  requirements 
into  the  computation  of  thrust  limit  on  the  growth 
distance.  Even  though  the  computation  of  thrust  limits 
on  the  growth  distance  can  be  performed  offline  for 
the  nominal  operating  conditions,  fast  computational 
procedures  are  beneficial  in  case  of  thruster  failures, 
degradations,  and  restrictions  on  thrust  directions  (e.g., 
caused  by  the  presence  of  other  spacecraft  nearby),  all 
of  which  can  lead  to  changing  constraints  on  thrust 
during  spacecraft  missions.  If  the  thrust  limits  are  pre¬ 
scribed  in  the  form  of  2-norm  bounds,  the  optimization 
problem  involved  in  computing  the  thrust  limit  on  the 
growth  distance  is  non-convex  (unlike  computing  the 
growth  distance  itself).  Consequently,  we  used  poly¬ 
hedral  norm  bounds  on  thrust,  that  lent  themselves  to 
explicit  solutions. 

In  this  paper,  we  address  the  problem  of  avoiding  a 
moving  debris.  Our  connectivity  graph  construction  uti- 
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lizes  the  sets,  centered  about  the  equilibria  in  the  virtual 
net,  that  are  guaranteed  to  contract  to  specified  smaller 
subsets  in  N  steps.  The  connectivity  of  two  equilibria 
in  the  net  ensures  that  over  N  steps  the  transition  to 
the  neighborhood  of  the  target  equilibrium  (i.e.,  smaller 
subset)  can  occur  without  the  debris  colliding  with  the 
spacecraft.  The  graph  vertices  are  formed  based  on 
equilibria  and  the  control  time  intervals  of  N  steps. 
A  graph  search  algorithm  is  employed  to  determine  a 
feasible  debris  avoidance  trajectory. 

The  related  literature  on  spacecraft  trajectory  opti¬ 
mization  with  obstacle/debris  avoidance  is  surveyed  in 

[I] .  Previous  research  addresses  topics  in  spacecraft 
trajectory  optimization  [6],  collision  avoidance  strate¬ 
gies  based  on  risk  assessment  [7],  the  use  of  artificial 
potential  functions  [8],  [9],  and  the  use  of  conventional 
and  mixed  integer  linear  programming  techniques  [10], 

[II] ,  [12],  [13]. 

The  paper  is  organized  as  follows.  In  Section  II 
we  discuss  the  linearized  models  used  to  represent 
spacecraft  relative  motion  dynamics.  In  Section  III 
we  review  our  approach  to  constructing  the  virtual 
net  based  on  a  set  of  forced  and  unforced  equilibria, 
and  using  this  virtual  net  for  debris  avoidance.  The 
procedure  to  compute  the  thrust  limit  on  the  growth 
distance  is  presented.  Simulation  results  are  reported 
in  Section  IV.  Finally,  concluding  remarks  are  made  in 
Section  V. 


dynamics  are  completely  controllable  from  Fy  input 
but  are  not  controllable  from  Fx  input.  The  out-of- 
plane  dynamics  are  controllable  from  Fz  input.  These 
dynamics  are  clearly  different  from  typical  ground 
robots. 

Assuming  a  sampling  period  of  AT  sec,  we  can 
convert  the  model  (1)  to  a  discrete-time  form 

X(t  +  1)  =  AX(t)  +  BU(t),  (2) 

where  X(t)  =  [x(t),  y(t),  z(t),  x(t),  y(t),z(t)]T 

is  the  state  vector  at  the  time  instant  t  G  Z+, 
U(t)  =  [Fx(t),  Fy(t ),  Fz(t)]T  is  the  control  vector 
of  thrust  forces  at  the  time  instant  t  G  Z+,  and  A  = 
exp(AcAT),  B  =  JQAT  exp(Ac(AT  —  r))drBc  are  the 
discretized  matrices  obtained  based  on  the  continuous¬ 
time  system  realization  (. AC,BC )  in  (1).  Alternatively, 
the  control  vector  U  can  represent  an  instantaneous 
change  in  the  velocity  of  the  spacecraft,  Av,  induced 
by  thrust,  with  an  appropriately  re-defined  T-matrix, 

"  0  0  0  " 

0  0  0 


0  1  0 
_  0  0  1  _ 

III.  Debris  Avoidance  Based  on  a  Virtual 
Net 


II.  Relative  Motion  Model 

The  spacecraft  relative  motion  model  represents  the 
spacecraft  dynamics  in  the  (non-inertial)  Hill’s  frame 
with  the  origin  at  a  target  location  on  a  nominal 
circular  orbit.  The  nonlinear  model  is  presented  in 
[1]  and  [2]  For  small  distances,  the  linearized  Hill- 
Clohessy- Wiltshire  (HCW)  equations  [15]  approximate 
the  relative  motion  of  the  spacecraft  on  a  circular  orbit 

aS  Fx 

x  —  3  n2x  —  2  ny  — 

mc 

Fy 

y  +  2  nx  =  — ,  (1) 

mc 

2  Fz 
Z  +  TiZ  —  - , 

mc 

where  Fx,Fy ,  Fz  are  components  of  the  external  force 
vector  (excluding  gravity)  acting  on  the  spacecraft.  The 
linearized  dynamics  account  for  differences  in  gravity 
between  the  spacecraft  and  nominal  orbital  location, 
and  for  relative  motion  effects.  The  spacecraft  relative 
motion  dynamics  in  the  orbital  plane  (x  and  y)  and 
in  the  out-of-orbital  plane  (z)  are  decoupled.  The  in¬ 
plane  dynamics  are  Lyapunov  unstable,  while  the  out- 
of-plane  dynamics  are  Lyapunov  stable.  The  in-plane 


Our  approach  to  debris  avoidance  is  based  on  utiliz¬ 
ing  constraint-admissible  positively  invariant  sets  [14], 
[16]  centered  around  the  spacecraft  forced  and  unforced 
relative  motion  equilibria.  A  finite  set  of  these  equilibria 
used  for  constructing  debris  avoidance  maneuvers  is 
referred  to  as  a  virtual  net.  Given  an  estimate  of  the 
debris  position,  we  build  a  connectivity  graph  that 
identifies  the  equilibria  in  the  virtual  net  between  which 
the  spacecraft  can  move,  with  guaranteed  collision-free 
motion  and  within  the  available  thrust  authority.  We 
then  employ  graph  search  to  determine  an  efficient  path 
between  the  equilibria  that  ensures  debris  avoidance. 

A.  Virtual  Net 

The  virtual  net  comprises  a  finite  set  of  equilibria, 
Xe(r),  corresponding  to  a  finite  set  of  prescribed  space¬ 
craft  positions  r  G  N  =  {rr,  7*2, . . . ,  rn}  C  i?3, 


Xe(rk)  = 


rkx 

rky 

rkz 

0 

0 

0 


k  —  1,  ■ 


,  n. 


(3) 
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whose  velocity  states  are  zero,  and  where  n  is  the 
number  of  equilibria  in  the  virtual  net.  See  Figure  1. 
We  assume  that  for  all  r  G  J\f,  the  corresponding  values 
of  control  necessary  to  support  the  specified  equilibria 
in  steady-state  satisfy  the  imposed  thrust  limits. 


Fig.  1:  The  virtual  net  for  debris  avoidance.  Dots 
correspond  to  positions  at  equilibria,  Xe(r),  on  a  virtual 
net.  The  ellipsoid  represents  the  debris  position  and 
uncertainty. 


B.  LQ  Controller  with  Gain  Switching 

A  conventional  Linear-Quadratic  (LQ)  feedback  is 
used  to  control  the  spacecraft  to  a  commanded  equilib¬ 
rium  in  (3), 

U  =  K[X  -  Xe(r))  +Tr  =  KX  +  H(K)r ,  (4) 


where 


T  = 


—3  n2mc  0  0 

0  0  0 
0  0  n2mc 


H(K )  =  r  -K 


h 

O3 


and  where  1 3  denotes  the  3  x  3  identity  matrix  and 
O3  denotes  the  3x3  zero  matrix.  This  LQ  controller 
provides  an  asymptotically  stable  closed-loop  system 
but  does  not  enforce  the  constraints. 

To  provide  greater  flexibility  in  handling  constraints, 
a  multimode  controller  architecture  is  employed  [16]. 
Specifically,  we  assume  that  a  finite  set  of  LQ  gains 
K  G  /C  =  {iTi,  •  •  •  ,  Km}  is  available  to  control  the 
spacecraft.  By  using  a  large  control  weight  in  the  LQ 
cost  functional,  motions  with  low  fuel  consumption 
yet  large  excursions  can  be  generated;  using  a  large 
control  weight  in  the  LQ  cost,  motions  with  short 
transition  time  can  be  generated  [17].  We  assume  that  a 
preference  ordering  has  been  defined  and  the  gains  are 
arranged  in  the  order  of  descending  preference,  from 
K\  being  the  highest  preference  gain  to  Km  being  the 
lowest  preference  gain. 


C.  Positively  Invariant  Sets 
The  ellipsoidal  set 

C(r,  K)  =  {X  €  R6  : 

\{X  -  Xe(r))TP(K)(X  -  Xe(r))  <  1}  c  R6,  (5) 

where  A(K)tPA(K)  -  P  <  0,  A(K)  =  (A  +  BK ), 
P  =  P(K)  >  0,  is  positively  invariant.  Positive 
invariance  implies  that  any  trajectory  of  the  closed-loop 
system  that  starts  in  C(r,K)  is  guaranteed  to  stay  in 
C(r,K)  as  long  as  the  same  LQ  gain  K  is  used  and  the 
same  set-point  command  r  is  maintained.  To  achieve 
the  positive  invariance,  the  matrix  P  can  be  chosen 
as  the  solution  of  the  discrete-time  Riccati  equation  or 
of  the  above  Lyapunov  equation  for  the  closed-loop 
asymptotically  stable  system.  We  note  that,  because 
the  system  is  linear,  the  positive  invariance  of  C(r,K) 
implies  the  positive  invariance  of  the  scaled  set 

C(r,  K,  p)  =  {X  €  R6  : 

Xe(r))T P{K)(X  -  Xe(r))  <  p2},  p  >  0. 

(6) 

Geometrically,  the  set  67(r,  iT,  p)  corresponds  to  an 
ellipsoid  scaled  by  the  value  of  p  and  centered  around 

Xe(r),  r  G  AT. 


D.  Debris  Representation 

We  use  a  set,  0(z,  Q ),  centered  around  the  position 
z  G  P3,  to  over-bound  the  current  position  of  the 
debris,  i.e., 


0(z,  Q)  =  {X  G  R6  :  (SX  -  z)tQ(SX  -  z)  <  1}, 

(7) 

where  Q  =  QT  >  0  and 
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1  0  0  0  0  0 
0  1  0  0  0  0 
0  0  1  0  0  0 


(8) 


The  set  0(z,  Q )  can  account  for  the  debris  and  space¬ 
craft  physical  sizes  and  also  for  the  uncertainties  in  the 
estimation  of  the  debris/spacecraft  position.  Note  that 
the  set  0(z,  Q )  has  an  ellipsoidal  shape  in  the  position 
directions  and  it  is  unbounded  in  the  velocity  directions. 
Ellipsoidal  sets  rather  than  polyhedral  sets  are  used 
to  over-bound  the  debris  since  ellipsoidal  bounds  are 
typically  produced  by  position  estimation  algorithms, 
such  as  the  Extended  Kalman  Filter  (EKF). 
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E.  Debris  Avoidance  Approach 

Consider  now  r*  G  A f,  representing  a  possible 
position  on  the  net  that  the  spacecraft  can  move  to  as 
a  part  of  the  debris  avoidance  maneuver.  Suppose  first 
that  the  debris  is  stationary  and  that  the  current  state 
of  the  spacecraft  is  X{to)  at  the  time  instant  to  G  Z+. 
If  there  exists  a  p  >  0  and  Kj  G  1C  such  that 

X(t0)  G  C{ruKj,p)  and  0(z,  Q)  n  C(ru  Kj,  p)  =  0, 

(9) 

the  spacecraft  can  move  to  the  position  r{  G  AT  by 
engaging  the  control  law  with  r(t)  =  rh  and  K(t )  = 
Kj,  t  >  to,  and  without  hitting  the  debris  confined 
to  0(z,  Q).  This  idea  underlies  our  approach  to  debris 
avoidance,  in  which  we  maintain  the  spacecraft  trajec¬ 
tories  within  the  tube  formed  by  positively  invariant 
sets  that  do  not  overlap  with  the  debris. 

To  avoid  a  non- stationary  debris,  in  [1],  [2]  we 
assume  that  its  path  can  be  covered  by  a  union  of  a 
finite  number  of  sets  of  ellipsoidal  shape, 

l — Tld 

D=  |J  0(zi,Qi),  (10) 

1  =  1 

where  the  center  of  the  Zth  set  is  denoted  by  zi  G  R 3 
and  the  Zth  set  shape  is  defined  by  Qi  =  Qj  >  0.  Then 
the  debris  path  avoidance  condition  is 

X(t0)  G  C(ri,  Kj,  p)  and 
0(zh  Qi)  n  C(ri,  Kj ,  p)  =  0,  l  =  1, . . . ,  nd. 

(11) 

The  same  approach,  with  larger  rid,  can  be  used  to 
handle  multiple  non- stationary  debris.  Note,  however, 
that  this  approach  is  conservative  as  it  does  not  account 
for  the  debris  progressions  along  their  paths  versus 
time. 

In  this  paper,  we  introduce  the  notion  of  time  into 
the  problem;  whereas  a  transition  between  r*  and  r3 
might  not  be  feasible  at  time  t\,  based  on  the  motion 
of  a  debris,  it  might  become  feasible  at  time  £2- 
To  accommodate  moving  debris,  we  introduce  a  set 
Cfc(r,  K,  p),  0  <  k  <  N,  defined  by  the  following 
relation, 

A(K)k(ck(r,K,p)-{Xe(r)}) 

V  x  7  (12) 

C(r,K,p)-{Xe(r)}j, 

Note  that  if  X(0)  G  Ck{r,  K,  p),  then  X(l)  G 
Ck^[r,K,p\  X(2)  G  Ck-2(r,K,p),-"  ,X(k)  G 
Co(r,K,p)  =  C{r,K,p).  The  set  Ck(r,K,p)  can  be 
much  larger  than  C(r,  iT,  p)\  any  states  in  Ck(r,  K ,  p) 
contract  to  C(r,K,p)  in  k  steps. 


F.  Connectivity  Graph  and  Graph  Search 

We  now  introduce  a  notion  of  connectivity  between 
two  vertices  of  the  virtual  net,  G  Af  and  r3  G  AT  at  a 
specified  time  to.  This  notion  is  based  on  the  fact  that 
the  time  to  transition  from  any  state  in  Cjv(r,  iT,  p)  to 
C(r,  iT,  p)  is  less  or  equal  than  N  steps.  Suppose  that 
the  debris  path  D(to  :  to  +  N  •  H)  has  been  predicted 
over  the  N  •  H  discrete  steps  from  the  time  instant  to, 
where 

t=tr 

D(tk  :tr)=  (J  0(z(t),Q(t)). 

t — tk 

The  node  T{  G  Af  is  connected  to  r{  G  Af  at  the  time 
instant  tk  =  to  +  kN  if  there  exists  K  G  1C  such  that 

D(tk:tk  +  N)nC(ri,K,p)  =  {®}.  (13) 

The  node  ri  G  Af  is  connected  to  node  r3  G  Af  at  time 
tk  if  there  exists  K  G  1C  such  that 

D(tk  :tk  +  N)n  CN(r3,K,  p)  =  {0}  (14) 

and 

C(ri,K,p)  cCN(rj,K,p).  (15) 

The  connectivity  implies  that  a  spacecraft  located 
close  to  an  equilibrium  corresponding  to  rz,  Xe(r^), 
can  transition  to  an  equilibrium  Xe(rj)  between  the 
time  instants  tk  and  tk  +  N  while  avoiding  collision 
with  the  debris.  We  note  that  if  is  connected  to  r3 
this  does  not  imply  that,  in  turn,  r3  is  connected  to  r% . 
We  also  note  that  connectivity  depends  on  the  existence 
of  an  appropriate  control  gain  from  the  set  of  gains  fC 
but  does  not  need  to  hold  for  all  gains.  Furthermore, 
since  connectivity  depends  on  the  predicted  motion  of 
the  debris,  connectivity/non-connectivity  can  depend  on 
time. 

The  on-line  motion  planning  with  debris  avoidance 
is  performed  according  to  the  following  procedure: 

Step  1:  Determine  the  debris  location,  shape  and 
predict  the  debris  path  D(to  :  to  +  N  •  H) 

Step  2:  Construct  graph  connectivity  matrices  cor¬ 
responding  to  tk,k  =  0, 1,  •  •  •  ,  H.  In  the 
graph  connectivity  matrix,  if  two  vertices,  r{ 
and  r3 ,  are  not  connected  at  tk,  the  corre¬ 
sponding  matrix  element  is  zero;  if  they  are 
connected  the  corresponding  matrix  element 
is  1.  In  parallel,  build  the  control  gain  selec¬ 
tivity  matrix,  which  identifies  the  index  of  the 
highest  preference  gain  K  for  which  ri  and 
Tj  are  connected.  This  gain  will  be  applied  if 
the  edge  connecting  and  r3  is  traversed. 

Step  3:  Perform  graph  search  to  determine  a  se¬ 
quence  r[tk\  G  Af  and  control  gains  K[k]  G 
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/C,  k  =  1  ,•••  ,/p,  such  that  r[t\\  satisfies 
the  initial  constraints,  r[lp\  satisfies  the  fi¬ 
nal  constraints,  and  the  path  length  lp  (or 
another  cost  function  such  as  the  expected 
fuel  consumption  or  expected  travel  time)  is 
minimized. 

Per  the  above  algorithm,  a  graph  search  is  utilized 
to  determine  the  minimum  number  of  equilibrium  hops 
around  a  debris  starting  at  to* 

Remark  1:  The  condition  (14)  is  conservative.  It  can 
be  replaced  by  a  less  conservative  condition, 

D(tk  :tk  +  m )  (~l  CN-m(rj,K,p)  =  {0}, 
m  =  0, 1,  •  •  •  ,  TV, 

at  a  price  of  more  demanding  computations. 

Remark  2:  The  condition  (14)  is  checked  com¬ 
putationally  using  the  fast  growth  distance  algorithm 
described  in  [1],  [2].  The  intersection  is  empty  if  Cjy 
can  be  grown  before  it  touches  D{tk  :  t&  +  N).  This 
fast  growth  distance  algorithm  is  essential  to  be  able  to 
rapidly  construct  the  connectivity  matrices. 

Remark  3:  In  our  simulations,  the  path  search  is 
performed  using  the  standard  Dijkstra’s  algorithm.  It  is 
applied  to  a  lifted  graph  with  vertices  being  the  pairs 
(?T  tk). 

IV.  Simulation  Results 

Simulations  are  now  provided  that  illustrate  our 
debris  avoidance  approach,  and  contrast  it  with  our 
prior  results  in  [1],  [2].  We  consider  a  nominal  circular 
orbit  of  850  km  and  discretize  the  HCW  equations  with 
a  sampling  period,  AT  =  30  sec.  We  construct  an 
approximately  2  km  cubed  virtual  net.  We  let  R  = 
2  x  10 7h  and  Q  =  diag( 100, 100, 100, 107, 107, 107), 
K  =  Ki. 

We  first  consider  the  case  of  a  non- stationary  de¬ 
bris  where  we  treat  its  motion  as  the  union  of  static 
debris  along  the  path  as  described  in  [2].  A  union 
of  ellipsoidal  sets  over-bounds  the  debris’  motion, 
U|=id  0(zi,Qi)9  where  Zi  are  generated  by  sampling 
the  relative  motion  of  the  debris  with  the  initial  condi¬ 
tion  [0  0.5  0  0  0.0006  0]T,  and  with  Qi  =  2OO/3,  i  = 
1  ...rid.  The  spacecraft’s  initial  condition  is  X(0)  = 
Xe(ro),  where  ro  =  [0  1  0]T.  The  target  equilibrium 
node  is  Xe(0).  We  impose  a  maximum  thrust  constraint 
of  10  N  in  each  axis.  Dijkstra’s  algorithm  is  used  to  find 
the  shortest  length  path  from  the  initial  node  to  the  final 
node.  The  simulation  results  in  Figure  2  demonstrate 
the  the  spacecraft  is  able  to  avoid  the  closed  debris 
path  by  hopping  under  it. 


We  now  repeat  this  simulations  using  the  approach 
described  in  this  paper  which  takes  into  account  the  de¬ 
bris’  motion  as  a  function  of  time.  Figure  2  shows  that 
the  graph  search  algorithm  is  able  to  find  a  path  which 
passes  through  the  debris’  path  but  avoids  collision  due 
to  the  debris’  location  elsewhere  at  the  specific  time 
instant  at  which  the  spacecraft  path  crosses  the  debris’ 
path. 


Fig.  2:  (a)  Debris  avoidance  path  for  a  non- stationary 
debris  using  the  union  method  described  in  [2].  The 
green  x  marks  the  initial  node.  The  blue  x  marks  the 
final  node.  The  red  ellipsoids  represent  the  debris  path. 
The  blue  line  is  the  path  the  spacecraft  takes  in  order 
to  avoid  the  debris.  The  blue  ellipsoids  represent  the 
maximally  grown  (see  [2])  invariant  sets,  C ,  along  the 
path,  (b)  The  time  history  of  thrust  magnitude. 


V.  Conclusion 

This  paper  described  an  approach  to  spacecraft  ma¬ 
neuver  planning  that  uses  contractive  constraint  admis¬ 
sible  sets  in  order  to  avoid  collisions  with  a  moving 
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(a)  Debris  Avoidance  Path. 


Fig.  3:  (a)  Debris  avoidance  path  for  a  non- stationary 
debris  using  the  contractive  set  approach  described  in 
this  paper.  The  green  x  marks  the  initial  node.  The  blue 
x  marks  the  final  node.  The  red  ellipsoids  represent  the 
debris  path.  The  blue  line  is  the  path  the  spacecraft 
takes  in  order  to  avoid  the  debris.  The  blue  ellipsoids 
represent  the  invariant  sets,  CV,  along  the  path,  (b)  The 
time  history  of  thrust  magnitude. 


debris.  Future  work  will  explore  cost  functions  other 
than  the  path  length  and  include  set-bounded  uncer¬ 
tainty  into  the  motion  planning. 
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